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ABSTRACT 

High signal-to-noise, representative spectra of star-forming galaxies at z ~ 2, obtained via stacking, 
reveal a high-velocity component underneath the narrow Ha and [Nil] emission lines. When modeled 
as a single Gaussian, this broad component has FWHM > 1500 km s _1 ; when modeled as broad 
wings on the Ha and [Nil] features, it has FWHM > 500 km s _1 . This feature is preferentially found 
in the more massive and more rapidly star-forming systems, which also tend to be older and larger 
galaxies. We interpret this emission as evidence of either powerful starburst-driven galactic winds or 
active supermassive black holes. If galactic winds are responsible for the broad emission, the observed 
luminosity and velocity of this gas imply mass outflow rates comparable to the star formation rate. 
On the other hand, if the broad line regions of active black holes account for the broad feature, the 
corresponding black holes masses are estimated to be an order of magnitude lower than those predicted 
by local scaling relations, suggesting a delayed assembly of supermassive black holes with respect to 
their host bulges. 

Subject headings: galaxies: high redshift 
active 



galaxies: evolution - galaxies: emission lines - galaxies: 



1. INTRODUCTION 

Observations of galaxies in the early Universe are a 
unique probe of matter assembly during its most active 
epoch; at z ~ 2, both the cosmic star formation rate and 
the luminous quasar space density are a t their peaks (e.g. 
iFan et all 120011: IChapman et al.l l2005h . Galaxies them- 
selves undergo corresponding growth during this time, 
with the total stellar mass density in galaxies increasing 
from ~ 15% to 50 —75% its current value between z ~ 3 
and z ~ 1 (e.g. iDickinson et alj l2003t iRudnick et al.1 
[200112001 . Constr aining the dynamical and baryonic 
processes driving this rapid evolution is therefore central 
to our understanding of galaxy formation and to inform- 
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ing cosmological simulations. 

At the relevant epochs, key spectral diagnostic features 
are redshifted into the near-infrared. In recent years, a 
number of surveys have therefore begun to systemati- 
cally probe high-r edshift population s wi th near-infrared 
spect r oscopy (e.g. lErb et alj l2006alfblfp: [Swinbank et al.l 
120041: iTakata et al.l l2006t iKriek et al.l l2008f) . Our k- 
cently completed SINS (Spectroscopic Imaging in the 
Near-infrared with SINFONI) survey has combined the 
resolving power of 8-10m class telescopes with high- 
resolution integral field spectrographs to study the de- 
tailed internal processes at wor k within massive, star- 
forming galaxies at z ~ 2 (iForster Schreiber et aLl 
20061 120091 iGenzel et al.ll2006l 120081 : iBouche et al.l l2007 



Shapiro i et al. I 20081: Cresci et al. 120091: see also re- 
lated work by IPuech et alj 120061: ISwinbank et all 120061 : 
Wright et alj|2007t lLaw et al.l 120071 12009ft . 

In this paper, we combine the spectra of 47 galaxies 
detected in Ha emission to study the average spectral 
properties of star- forming galaxies at z ~ 2 (^). In par- 
ticular, we report the discovery of broad emission lines in 
these galaxies ([J3]) and interpret this high- velocity warm 
gas as arising either in large-scale galactic winds driven 
by the high star formation rates (SFR) in these galaxies 
or in the broad-line regions (BLR) surrounding active 
galactic nuclei (AGN). We explore the implications of 
both scenarios in g] and conclude in fj5] 
Throughout this paper, we assume a A-dominated 



cosmology with Hq — 70 km s Mpc 



0, 



0.3, 



and J1a = 0.7. For this cosmology, 1" corresponds to 
w 8.2 kpc at z = 2.2. 

2. DATA AND ANALYSIS 

In the context of the SINS program, 80 z =1—3 
systems were observed in emission lines in the 
infrared (rest-frame optical) with VLT/SINFONI 
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(|Eisenhauer et al.ll2003l ; iBonnet et al.ll2004D for an aver- 
age of 3.5 hours per band and pixel scale on each target 
(jForster Schreiber et al.l I2009D . These galaxies were 
largely (62/ 80) taken f r om the (rest-frame) UV-selected 
samples of lErb et aU ( 2006b. c ) and the (res t- frame ) 
optically-sele c ted sa m ples of [Abraham et al.l (f2 004). 
DaddGTaD (I2004bl) . iKong et all (|2006ft . iLillv et all 
( 20071 ). and Kurk et al. (in prep). From these 
magnitude- and color-defined samples, suitable SINS 
targets were culled, with the main selection criteria 
being a combination of target visibility during the 
observing runs, night sky line avoidance in the emission 
lines of interest, and an estimated integrated emission 
line flux > 5 x 10~ 17 erg s" 1 cm~ 2 , such that high 
quality data could be obtained in reasonable integration 
times. Of the 62 rest-frame UV/optically-selected 
galaxies chosen in this manner, 52 were well detected in 
H a in our SINFONI obse r vatio ns . 

IForster Schreiber et a l. (2009) discuss the selection of 
this sample in detail and show that the SINS galaxies 
are representative of the z ~ 2 star-forming galaxy pop- 
ulation, with some bias towards the more rapidly star- 
forming (and therefore more luminous in Ha emission) 
systems. They also note that the SINS sample (and some 
of the parent samples) selects against known AGN and 
quasars, in the interest of studying the dynamic and evo- 
lutionary state of z ~ 2 star- forming galaxies. However, 
a small number of previously known AGN were in fact ob- 
served in the SINS program; in the UV/optically-selected 
part of the sample, there are 5 such systems, as origi- 
nally identified with the UV/optical spectroscopy of the 
parent surveys. For further details about the SINS sam- 
ple, obser vations, and data reduction, we refer interested 
readers to IForster Schreiber et alJ (|2009h . 

Here, we analyze our SINS observations of the galax- 
ies that were UV/optically- identified and that are well- 
detected in Ha in individual spatial elements; this pop- 
ulation (totalling 47 of 52 sources detected in Ha, in- 
cluding 4/5 of the previously known AGN) comprises 
the majority of the SINS sample. To study the average 
properties of these z ~ 2 star-forming galaxies, we gen- 
erated stacked, high-signal-to-noise (S/N) spectra repre- 
sentative of the population as a whole. 

From our integral field data, we first created a 
spatially-integrated one-dimensional spectrum for each 
galaxy by shifting each spectrum within the galaxy dat- 
acube by its measured (Ha) velocity and then collaps- 
ing the datacube into a single spatially-integrated spec- 
trum. In this manner, the spatially- integrated spec- 
trum contains no systematic velocity broadening (e.g. 
by large-scale rotation) on scales larger than the PSF 
(FWHM ~ 4 kpc). Testing confirmed that this approach 
did not affect the properties of the broad Ha component 
(see below) but did improve the signal-to- noise (S/N) of 
the detection. This technique has the additional benefit 
of randomizing OH atmospheric emission lines in the Ha 
rest-frame and therefore effectively eliminating residuals 
from the OH line removal. The remaining residuals from 
this process were inspected and masked out by hand. 

The 47 spatially-integrated spectra were then com- 
bined into a single spectrum (with equivalent integra- 
tion time of 195 hours) by interpolating all spectra onto 
a common wavelength axis, converting their measured 
fluxes to luminosities using their luminosity distances, 



weighting each spectrum by the S/N of the Ha emis- 
sion line, and averaging the resulting spectra. During 
this process, we do not correct for extinction (but see 
£j4j) . Typical extinctions o ur sample galaxies have been 
measured t o be Ay ~ 1 jF orster Schreiber e t al.l [2009L 
see also e.g. iDaddi et al.ll2004at lErb et al.ll2006d) . which 
translates to an underestimation of our Ha luminosities 
by at most a factor of ~ 2. 

The average spectrum for the SINS z ~ 2 star-forming 
galaxies is presented in the top panel of Figure [TJ We 
also created average spectra of subsets of the SINS galaxy 
sample, in order to test the dependence of spectral prop- 
erties on other known galaxy properties. 

Our average spectra reveal a broad emission compo- 
nent underneath the bright narrow lines. We quantify 
this feature in each average spectrum by simultaneously 
fitting a combination of a constant continuum offset, nar- 
row lines (Ha, [Nil], [SII]) of identical kinematics (veloc- 
ity and velocity dispersion), and a single broad compo- 
nent, whose kinematics are allowed to vary. All lines are 
assumed to be well-described by a single Gaussian; the 
validity of this assumption is confirmed by a reduced x 2 
of close to unity (x 2 do f ~ 0.9-1.9) for all fits (see TableQJ. 
Our data can also be well fit by fitting a combination of a 
constant continuum offset, narrow lines with shared kine- 
matics, and broad forbidden and permitted lines with 
shared kinematics, with the [Nil] /Ha ratio identical in 
the narrow and broad component. All lines are assumed 
to be well-described by this combination of two Gaus- 
sians (dotted green line in Figure [TJ Xdof = 1-4), which 
has the same number of free parameters as the previous 
fit. In the limit of the S/N of our data, we cannot add 
additional free parameters to the fits, nor can we identify 
a preferred model. For clarity, we refer throughout to the 
former (single broad Gaussian under the Ha+[NII] com- 
plex) as "broad lines" and the latter (double Gaussians 
for both Ha and [Nil]) as "broad wings." For simplic- 
ity and for comparison with the literature, we primarily 
quantify the observed high-velocity feature with a sin- 
gle broad Ha line in fj3j However, we also discuss the 
implications of the broad wings scenario on the derived 
properties of the broad emission (fj3j> and on the inter- 
pretation of this emission (SJU). 

From these fits to the average spectra, we measure the 
fractional contribution of broad emission to the overall 
emission line flux, the kinematics of the broad compo- 
nent, and the line ratios of the narrow line components. 
The significance of these measurements arc quantified 
in two ways. First, for each average spectrum, we re- 
created 100 spectra by randomly sampling (with replace- 
ment) and combining the individual contributing galaxy 
spectra. The properties of the broad component were 
measured in each of the resulting 100 average spectra, 
yielding confidence intervals for all derived quantities. 
Second, the probability of false positive detections F 'false 
was tested by creating 1000 simulated spectra with the 
narrow line and observational properties (noise, spec- 
tral resolution) characteristic of each average spectrum. 
Comparing the derived broad components in these spec- 
tra to those in the real SINS spectra, we estimate the rate 
of spurious detections of broad components equal to or 
more prominent (in luminosity and FWHM) than those 
in the actual data. The results of this analysis for the 
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Fig. 1. — Left: Average spectrum of the 47 SINS galaxies (top panel), followed by the average spectrum of all non-active systems and 
that of all systems previously known to host AGN (lower two panels). Right: Zoomed view of the He? and [No] region, on a velocity scale, 
with the best-fitting combination of a constant continuum, narrow lines, and a broad component overplotted (red). The individual line 
components are also plotted (blue). For reference, also overplotted are the best fits derived by fitting both Ho and [Nn] with only narrow 
lines (dashed magenta) and with a narrow and a broad component (dotted green). 



different average spectra are presented in the following 
section. 

3. RESULTS 

We find that the average SINS galaxy spectrum in- 
cludes a significant amount of broad emission (top panel 
of Figure [l}, with Xdof = 4-9 for a fit with only narrow 
lines (dashed magenta line in top right panel of Figure [1} 
and Xdof — 1-8 an d 1-4 for fits including a broad Ha line 
and broad Ha and [Nil] wings, respectively (red and dot- 
ted green lines in same panel; P 'false = 2% and 1%). To 
ensure that this signature is not the result of the 4 known 
AGN included in this sample, we also create stacked spec- 
tra of the AGN and the rest of the sample. While the 
broad emission from the AGN host systems alone is quite 
substantial, a comparison of the average SINS spectrum 
and the average non-AGN spectrum illustrates that the 
broad emission in the average SINS spectrum is not dom- 
inated by that coming from the 4 AGN. These results are 
summarized in Table [T] 

We test the dependence of the presence of broad 



emission on galaxy properties by dividing the sam- 
ple into three stellar mass bins, using the results 
from the spectral ener g y dist ribution (SED) fitting of 
iForster Schreiber et all (20091 ) for the 45 of our 47 tar- 
gets for which sufficient broadband data exist. The aver- 
age spectra for these three bins show an increasing pres- 
ence of a broad component with stellar mass (Figure [51 
Pfaise — 11%, 5%, 3% respectively). However, the spec- 
trum of the highest mass bin is significantly affected by 
the contribution of the 4 previously known AGN in our 
sample, which all fall into this bin. The nuclear emission 
in these systems can bias the results of SED fitting to- 
wards larger masses, so we confirm their high masses with 
the dy namical mass measurements made bv lCresci et al.l 
( 2009); in all cases, the dynamical masses of these galax- 
ies are consistent with the stellar masses used here and 
remain among the highest in the SINS sample. Never- 
theless, we confirm that a (weaker) broad component is 
also present in the non-active galaxies in this bin (green 
line in Figure [2]), with P 'false = 9%. 
Several other key properties of galaxies, including SFR, 
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Fig. 2. — Left: Average spectrum of each mass bin. Right: Zoomed view of the Ha and [Nn] region, on a velocity scale, with the 
best-fitting combination of a constant continuum, narrow lines, and a broad component overplotted (red). The individual line components 
are also plotted (blue). The best fit obtained with only narrow emission lines is plotted for comparison (dashed magenta). In the high mass 
bin, the best-fit broad line component to the average spectrum of the 3 non-active high mass galaxies is also overplotted (dotted green). 



size, stellar age, and metallicity, have well established 
correlations with stellar mass in high-z galaxies (e.g . 
iNoeske et aI1l2007t iTruiillo et al.ll2006b lErb et al.ll2006rj ) . 
Both the SFR-M* relation and the mass-metallicity rela- 
tion are apparent in Figure via the increasing narrow 
Ha luminosity and increasing [Nil] /Ha ratio with stellar 
mass, respectively. Since the [Nil] /Ha ratio remains well 
below levels expected of shock heating or AGN activity, 
this latter is most likely tracing variations in metallicity 
(see below, as well as Buschkamp et al. in prep). To these 
established mass-dependent properties in high-z galax- 
ies, we now add the presence and strength of a broad 
component. 

With the spatially resolved data, we can also compare 
the integrated spectra from the central (R < 3 kpc) re- 
gions of high-z galaxies to those from extended (i?, = 3 — 
15 kpc) regions, in order to determine what regions in 
these systems are generating broad emission. For this 
analysis, we use galaxies from the intermediate and high 
mass bins of Figure in which the intensity distribution 
of the stellar continuum defines a clear center of the sys- 
tem (totalling 6 systems). The average spectra of the 



central and extended regions of these systems are shown 
in Figure [3J normalized to the spatial area over which 
the spectra were extracted. 

In these spectra, a broad component is preferred by 
the best-fitting models; however, the significance of this 
result is low. The detection of broad emission in galaxy 
centers at z ~ 2 (with P/ Q / se = 8%) is somewhat more 
robust than that in the extended regions (P false = 14%), 
in which the best fit broad component is very shallow. If 
real, the broad feature in the extended regions accounts 
for a comparable fraction of the total Ha luminosity to 
that in the central regions (Table [l}. Tests of simulated 
galaxies indicate that such a broad line in the extended 
regions cannot be reproduced by a nuclear point-source 
of broad emission (i.e. AGN) broadened by the PSF. 

In the spectra shown in Figure [3J we note that the 
average [Nil] /Ha ratio is comparable in the central and 
extended regions, despite the difference in broad emission 
component, in support of the above interpretation of the 
[Nil] /Ha feature as primarily reflecting the metallicity of 
these systems. However, this should not be interpreted 
as a lack of metallicity gradient in these galaxies, since 
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Fig. 3. — Average spectra, in luminosity surface density, of the central (top panel) and extended (bottom panel) regions of well-resolved 
SINS galaxies. Panels and colors are as in Figures [T] and f2] 



TABLE 1 
Results from Line Fitting 



Subsamplc 


Figure 
Reference 


L H OL,broad 

(10 41 erg s" 1 ) 


f a 
j broad 


v.// 6 
(km s ) 


FWHM„ a ,i, roarl 
(km s — 1 ) 






X-do f , narrow 


All 


1, top 


4 0+ 11 


n 90+O.O4 


18+™ 


1632 +44b 


0.02 


1.8 


4.9 


All, broad Ho and [No] d 


1, top 


5 9+ 10 


42+ - 08 
u '* z -0.08 


si 3 4 


556 +2 s f 


0.01 


1.4 


1.9 


Non-AGN 


1, center 


, Q+0.6 


n 99+O.O3 
u - z -0.02 


1Q+82 
la -42 


1329+^J; 


0.03 


1.9 


3.7 


AGN 


1, bottom 


o O+0.7 
8 - z -2.8 


u - JD -0.09 


42 +94 
^ z -29 


2863 + «| 


0.007 


0.9 


3.0 


M* < 2 x 10 10 M 


2, top 


n 9+0.4 
^•^-0.6 


21+ 03 
u ' zl -0.02 


9 fi +86 


1051+/J 


0.11 


0.9 


1.3 


M* = 2 - 7 x 10 10 M Q 


2, center 




n 99+O.O6 
u - z -0.02 


in+71 
ii -70 




0.05 


1.4 


2.0 


M* > 7 x 10 10 M 


2, bottom 




3l+ () - 04 


9Q+82 


2183+|| 


0.03 


1.3 


3.6 


Center 


3, top 


r 1+2.8 e 

o.i_ 1 7 


24+ - 13 
u ' z ^-0.04 


O+80 
a -27 


1564+11 


0.08 
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1.9 


Extended 


3, bottom 


i n + 0.6 e 
i - y -0.5 


n 99+O.O6 
u - z -0.02 


10 +69 
iU -39 


1508_ 244 


0.14 


0.9 


1.2 



a Ratio of Ha luminosity in broad component to total Ha (narrow+broad) luminosity. 

6 Velocity offset of broad Ha feature from narrow Ha. 

c Reduced \ 2 for a fit assuming only narrow emission lines. 

d Most components of the fit (f //, FWHM^ Q 5 T . oad ), and fbroad) are constrained in the fit to be identical for the permitted and forbidden 
lines. The luminosity (L Ha broad) ls quoted only for the broad component of the Ha line. 
e Units are 10 39 erg s — 1 kpc — 2 . 



the spectra shown here are averaged over a number of 
galaxies, in each of which the extended region is spatially 
integrated over a large range in radii. Detailed studies of 
metallicity gradients within and between individual SINS 
systems will be presented in Buschkamp et al. (in prep). 

4. DISCUSSION 

The low luminosity broad emission seen in our SINS 
galaxies can be interpreted in one of two ways: either 
as evidence of large-scale galactic outflows, presumably 
driven by the galaxies' very high star formation rates, or 



as a tracers of the broad-line regions surrounding AGN 
in these early galaxies. The current data do not allow 
us to robustly distinguish between these two scenarios 
empirically; in the following, we therefore examine both 
in detail. 

4.1. Broad Emission from Starburst- Driven Winds 

One possible explanation of the broad emission is a 
starburst-driven wind. This scenario is in keeping with 
the positive correlation between broad emission and stel- 
lar mass (and therefore star formation rate) seen in Fig- 
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ure [2] and with the possible broad emission from non- 
nuclear regions seen in Figure El Moreover, starburst- 
driven winds are expected to be ubiquitous in the rapidly 
star-forming populations common at high-redshift (e.g . 



Pettini et a l. 2001; Sha plev et al. ! 2003; S mail et alJ2003t 



Weiner et al.ll2009i ) and are probably expelling mass from 



their host galaxies at rates comparable to the star for- 
mation rate (below; see also lMartinlll999t IPettini et al.l 
[20001: lMartinll200l lErbl[2008l ; 1Weiner et al.lT2009D . 

At low redshift, star- forming systems are known to 
drive galactic winds with observable signatures in the 
wings of the permitted an d forbidden emission li nes. In 
dwarf starburst galaxies, IWestmoquette et"aTI (|2007ah 
find that the broad wings of the emission lines can be 
modeled as a second Gaussian component with FWHM 
< 300 km s . In contrast, observations of the more mas- 
sive and more rapidly star-forming IR-luminous galaxy 
population reveal higher F WHM in the broad wings 
(300-800 km s'^ lArribas et al.ll2001h . whose line widths 
and f broad are comparable to those o bserved in the SINS 
galaxies (550 km s" 1 ; see Tabl e [Q lArmus et ail 119891 
Il990t ILehnert fc HeckmanI fl996h . 

Indeed, the relationship between wind speed and SFR 
(and the equivalent properties, galaxy mass and B- 
band magnitude) has b een d emonstrated by IRupke et al.l 
(|2005h . IMartinI (|2005h . an d iTremonti et alJ (I20Q2D, re - 
spe ctively. In particu lar, ILehnert fc HeckmanI (|1996l ) 
and lRupke et a l. ( 2005) showed that systems with higher 
SFR have faster winds, whose velocities exceed ~ 1000 
km s _1 in ionized gas tracers. We recover similar trends 
and velocities with the three mass bins of the SINS data, 
as shown in Figure |U Note, however, that this Figure 
should not be directly compared with those in the works 
listed above, in which the wind speeds are probed with 
interstellar absorption features, whose velocities are typ- 
ically lower than th o se of ionized emiss ion line gas (e.g. 
IVeilleux et all 120051 : IRupke et all l2005h . Nevertheless, 
the similarity between the properties of local winds and 
those of the broad wings seen in the SINS galaxies makes 
it plausible that this emission is due to the galactic winds 
that certainly exist in these systems. 

These broad wings can be produced via several mech- 
anisms, including turbulent mixing of the hot wind 
fluid and cool gas clumps at the base of the wind 
(|Westmoauette et al.l l2007al icT) or within t he large-scale 
bubb l e along the galaxy minor axis (e.g. lArmus et al.l 
119901: lArribas et a"l1l2001h . In the idealized model of 
a wind plowing into a s pheri cally symmetric ambient 
medium, iHeckman et al.l (|1993l ) provide analytic expres- 
sions for the velocity and Ha luminosity of the shocked 
gas as functions of the bolometric luminosity of the star- 
formation event, the duration of the event, and the ISM 
density. Converting the SFR (~ 100 M Q yr^ 1 ) of the 
SINS galaxies to a bolometric luminosity and accounting 
for the correspondingly dense r ISM than in loca l spiral 
galaxies (by a factor of 10— 30: lBouche etPa l. 2007), these 
relations predict a Ha luminosity of 0.7—10 x 10 42 erg 
s _1 generated from gas moving at ~ 70—300 km s . 
These properties are broadly consistent with the high- 
velocity wings observed in the SINS galaxies. 

However, our spectra can also be fit by a broad Ha line, 
and we therefore explore whether such a feature could 
likewise be generated in galactic winds. For simplicity, 
we assume that these winds are powered by supernova 
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Fig. 4. — FWHM of the broad component in the SINS galaxies 
in the three stellar mass bins defined in Figure [2] as modeled by a 
single broad Ho line (filled circles) and by broad wings on both the 
Ha and [Nil] lines (open circles). In the latter model, the low mass 
bin is not well constrained due to the low S/N of the [Nil] feature 
(see Figure [2}, resulting in a fit that strongly prefers a broad line 
for only the high S/N (Ha) line; we therefore omit this bin from 
the plot. 



remnants (SNR), in local examples of which broad Ha 
(FWHM = 500-8000 km s" 1 ) is observed throughout 
the Sedov- Taylor expansion phase, the result of charge 
exchange of the electrons fr om slow neutral atoms to the 
fast post-shock pro t ons (e.g.lChevalier fc Ravmondlll978l : 
iSmith et alJ fl99lL iHeng fc Sunvaevl l2008r ). Although 
such broad Ha emission is not observed on the galactic 
scale in local starbursting systems, we note that none of 
these systems are appropriate analogs to the star forma- 
tion mode that dominates at z ~ 2. Only local ULIRGs 
have SFR comparable to that in the SINS galaxies, and 
the ionized gas emission in ULI RGs suffers signific antly 
more extinction (Ay = 5— 1000: lGenzel et al.l ll998) than 
in SINS galaxies (Ay ~ 1). We therefore briefly examine 
the possibility that galactic winds at z ~ 2 emit broad 
Ha lines via a superposition of SNR. 

A simple test of this scenario is whether there are 
a sufficient number of SNR in the SINS galaxies to 
drive the observed broad Ha luminosity. Locally, 
SNR are observed to have broad Ha luminosities of 
10 30 — 10 34 erg s -1 , generated by change exchange as the 
shock encounters the ISM and therefore proportional to 
the square of the gas density. In the SINS galaxies, the 
high SFR (~ 100 M Q yr" 1 ) and dense ISM (10-30 times 
denser than in local spiral galaxies) imply an increase in 
broad Ha luminosity per SNR by a factor of 100—1000, to 
10 32 -10 37 erg s" 1 . The SFR in the SINS galaxies yields 
~ 1 supernova explosion per year, each of which will have 
expected Sedov- Taylor lifetimes of ~ 10 4 yr in the dense 
ISM, resulting in roughly 10 4 SNR radiating broad Ha 
at any given time in the average SINS galaxy, for a total 
broad Ha luminosity of 10 36 — 10 41 erg s _1 . This number 
is lower than the observed broad Ha luminosity in the 
SINS galaxies by a factor of a few (see Table QJ; however, 
real galactic winds penetrate much further through their 
host galaxies' ISM than do the sum of individual SNR, 
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implying that additional broad Ha emission is expected 
from the interaction of the large-scale winds with the 
ambient medium. This effect brings the predictions into 
even closer agreement with the observations and thus 
makes this mechanism an energetically feasible explana- 
tion for the luminosity and FWHM of the broad emission 
in the SINS galaxies. 

Regardless of whether the broad emission is a broad Ha 
line or broad wings on all emission lines, if the emission 
is in fact due to starburst-driven winds, it is worth exam- 
ining what the fate of this high- velocity gas may be. The 
average SINS galaxy (top panel of Figure [1]) has cither 
broad wings in Ha and [Nil] with velocity dispersions of 
~ 250 km s _1 or a broad Ha line with velocity disper- 
sion ~ 700 km s _1 . We can compare these values directly 
with the escape velocity v esc w 450 km s _1 for a typical 
SINS galaxy ((M dyn ) = 8 x 10 10 M , (R 1/2 ) = 3.4 kpc; 
iForster Schreiber et al.|[2009f ). A significant fraction (7% 
or 55%, respectively) of the high-velocity gas has veloc- 
ity exceeding the host galaxy's escape velocity; assuming 
that this gas is distributed throughout the star-forming 
disk, a non-negligible amount of it should be expected 
to escape the galaxy. Furthermore, if the surrounding 
dark matter halo has a flat rotation curve, this gas would 
also be expected to escape the halo into the intergalactic 
medium. 

We can then estimate the mass outflow rate that 
would correspond to such superwinds. The escaping Ha- 
emitting gas is ~ 3—17% of the total Ha-emitting gas 
in the SINS galaxies (i.e. 30% of the emission is broad 
and 55% of the broad emission escapes). Assuming that 
the fraction of gas in the ionized phase is roughly the 
same in the galaxies' star-forming disks and in the out- 
flows ( compare ~ 1% in th e Milky Way to ~ 0.1 — 1% in 
winds: IVeilleux et al.ll2005h . this implies that ~ 3—17% 
of the galaxies' gas reservoirs are being expelled by the 
observed star-forming event. The dynamical times asso- 
ciated with these outflows can be approximated as the 
radius of the star formation event divided by the veloc- 
ity of the flow; for the SINS galaxies, this yields dynam- 
ical times of ~ 10 Myr. The average SINS galaxy has 
a dynamical mass of 8 x 10 10 and a gas f raction 
of 0.2-0.4 (e.g. lErb et all 120063 : iBouche etall l200l . 
yielding an expected outflow rate of 50—5 yr" 1 . 
This is consistent with the results of lErbl |20Q8), who 
argue that the observed z ~ 2 mass-metallicity relation- 
ship requires outflow rates slightly larger than the SFR 
(SFR = 1-8 00 M m yr" 1 with median 72 M n yr" 1 in the 
SINS sample: [FoTster Schreiber et alj|2009h . If the broad 
emission we observe is due to galactic winds at z ~ 2, we 
thus find mass outflow rates consistent with the observed 
mctallicity evolution of these galaxies. 

4.2. Broad Emission from Active SMBHs 

Another possible, and common, interpretation of broad 
emission lines are as signatures of nuclear activity. For 
our galaxies, this interpretation is supported by the 
strength of the broad Ha emission increasing with stellar 
mass (and therefore bulge mass and possibly black hole 
mass) seen in Figure^ by the less luminous broad emis- 
sion in the non-active systems (the driving mechanism 
being "turned off" or obscured; Figures Q] and [2]) , and 
by the more significant detection of broad emission in 
the centers of galaxies than in their non-nuclear regions. 



Moreover, studies of z ~ 2 star-forming galaxies in in- 
frared and X-ray emission suggest that activ e nuclei may 
be a c ommon feature of this population dDaddi et al.l 
I2007al) . 

Early work in the local Universe has shown that broad 
Ha is by definition omnipresent in Type 1 AGN, in which 
the broad-line regions around the nuclei are unobscured. 
Moreover, the kinematics and sizes of these regions can 
be used to infer "virial" supermassive black hole (SMBH) 
masses, via calibrated relations between the observed 
AGN continuum l uminosities at 510 A (L5100) and the 
sizes of the BLR (jKaspi et al.ll2000D . which can then be 
used in combination with the width of the broad emis - 
sion lines to e s timat e SMBH masses (jVesterga ard 2002). 
iGreene fc Hoi (|2005l) have additionally related -L5100 to 
the luminosity of the broad emission lines, making it 
possible to measure virial black hole masses using only a 
single broad line, 



M BH = (2.0±g;|) x 10 6 

T \ 0.55±0.02 



10 42 erg s" 1 



FWHMhq 
10 3 km s- 1 



2.06±0.06 



M« 



(1) 

Such virial black hole mass estimates have been ver- 
ified against stellar and gas kinematic determinations 
of Mbh and found to be accurate to within a fac- 
tor^ ~3 (|Onken et ai] 120071: iHicks fc Malkanl 120081: 
iNetzerl [20091 but see also iMarconi et al.l I2008TT As a 
result, they have been utilized in the high- z Universe to 
probe the masses of SMBHs in quasars (e.g. Willott et al.l 



I2003t iMcLure fe DunlorJl2004 : IVestergaardll2004D and in 
the s ub-mm-bright galaxy population ( Alexander et al.l 
l2008h . 

Additionally, the broad-line Ha luminosity can also 
be used to probe the accretion rate of SMBHs. As- 
suming the AGN continuum luminosity L5100 is roughly 
1/10 the AGN bolometric luminosity iboi, the accre- 
tion rate can be es timated by Mbh = £boi/(?7c 2 ) 
(|La Mura et alJl2~007| ). Using a typical value of r\ ~ 0.1 
and the calibrated relationship between Lhq and L5100 
(| Greene fe Hoi I2005T ) . we can then estimate Mbh- Al- 
though approximate, this estimate of Mbh and the esti- 
mate of Mbh derived using equation Q] nevertheless allow 
us to study both the putative black hole masses in our 
SINS galaxies and their accretion rates using our mea- 
surements of the broad Ha features. 

With this procedure, we derive for the average SINS 
spectrum (top row in Figure [T|) a black hole mass of 



= 4+1 x 10 e M 
Furthermore 
9 x 10 5 

hi 
-0.5 



with an Eddington ratio of 
we infer black hole masses of 



+2 



10 6 M m , and 



M BH 
- 0.2 

M BH < 9 x 10 5 M Q , M BH = 2_ 2 a ±u m , 
Mbh = l^o 5 x 10' M Q in each of the three stellar mass 
bins, respectively. In the latter two bins, the Eddington 
ratios are estimated to be 0.3 and 0.4, respectively. 

We note that the Ha luminosities used in these black 
hole mass estimates have not been corrected for extinc- 
tion, and any such correction would increase the derived 
masses. However, broad emission lines associated with 
Type 1-1.5 AGN in local galaxies are consistent with no 
additi onal extinction along the line-of-sight to the BLRs 
(e.g. iRhee fe Larkh] 12000: lAlonso-Herrero et alj 120031: 
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Fig. 5. — Left: Estimated black hole masses (derived with v irial estimates from the broad He* emission line) against average stellar 
mass (derived with SED modeling; Forstcr Schreibcr ct al. 2009|) in each mass bin. Cent er. Relation between central mass concentration 
(dynamical mass within 3 kpc) and stellar mass for the five galaxies modeled in detail by Gcnzcl et al. (2008), along with the best-fitting 
third-order polynomial. Right: Estimated black hole masses against estimated bulge masses (circles; Mt, u i ge estimated from stellar masses 
shown in left panel and polynomial shown in cen ter panel) for the three mass bins shown in Figurc[2] Ovcrplotted (square) is the average 
BH mass measured in z ~ 2 sub-mm galaxies by Alexander et al. (2008); for c omparison w i th ou r central mass concentration (dynamical 
mass within 3 kpc), we adopt the stellar mass within ~ 4 kpc measured by Alexander et al. (2008) as the "bulge" mass, with the error bar 
extending to include all "bulge" masses estimated by those authors. The MBH-Mt, u i gc relation for the local Universe and its scatter are 
indicated with the solid line and shading (Haring & Rix 2004). 



I Greene fe" Ho 2005). Likewise, at z ~ 2, 1 Alexander et aTl 
( 2008f ) have used the Ha/H/3 broad- line Balmer decre- 
ment to measure only small amounts of nuclear extinc- 
tion (Ay ~ 1.2) in their dust-rich sub-mm population. 
These authors postulate that the plentiful dust in the 
sub-mm galaxies obscures regions of star-formation and 
not the BLRs. In the SINS galaxies, the average extinc- 
tion in the star- forming regions is An a ~ 0.8; if the BLR 
is similarly obscured, correcting for this effect would yield 
an increase of at most a factor of 2 in broad Ha lumi- 
nosity and, from equation [1] a factor ~ 1.5 in derived 
black hole mass. This suggests that it is unlikely that 
the black hole masses measured here suffer significantly 
from extinction. 

It is then naturally interesting to compare the esti- 
mated black hole masses with large-scale galaxy proper- 
ties. In the left panel of Figure [5l we find the expected 
trend that black holes of increasing mass are found in 
galaxies of increasing stellar mass, with black-hole-to- 
stellar mass ratios of ~ 7 x 10~ 5 . For such black holes to 
be consistent with local scaling relations between black 
holes and bulges (MeH/Mbuigc ~ 10~ 3 ), the bulge-to- 
total ratio B/T of these systems would need to be quite 
small (< 0.07). This is in marked contrast to results from 
detailed dynamical modeling of the five SINS galaxies 
with the highest-resolution observations (jGenzel et al.l 
2008), which yield B/T = 0.15-0.4 (center panel of Fig- 
ure [5j) . Parameterizing these results with a simple poly- 
nomial, we estimate "bulge" masses (~ dynamical mass 
within 3 kpc) for the three mass bins, albeit with large 
uncertainty. The resulting B/T for these mass bins are 
< 0.1, < 0.4, and 0.4, respectively. These bulge masses 
are then plotted against our measured black hole masses 
in the right panel of Figure [5] 

Comparing to the loc al MeH-Mbuige relation measured 
bv lHaring fc Rixl lj200l . we find that our z ~ 2 galaxies 
lie significantly below the local relation, implying that 
black holes in the star-forming galaxies in the early Uni- 
verse may have la gged significantly b ehind their host 
bulges in assembly. iGenzel et al.1 (|2008l ) have shown that 
z ~ 2 is, for many of these systems, an era of bulge forma- 
tion via smooth but rapid secular processes, during which 



massive bulges are assembled from large star-forming 
clumps (M ~ 10 s - 10^ M Q ) on timescales of < 1 Gyr. 
lElmegreen et al.l (|2008f ) have demonstrated in simula- 
tions that, assuming each such clump contains a black 
hole of 10 -3 of its total mass, these black holes would 
migrate to the galaxy center with their host clumps and 
form central SMBHs that are somewhat under-massive 
for the resulting bulges. The location of our z ~ 2 sys- 
tems significantly below the local relation likewise sug- 
gests that the bulges in these galaxies form first, through 
rapid secular processes, with the assembly of the central 
SMBHs following later. 

The timescale for this final SMBH growth can be esti- 
mated for the high mass bin, in which the galaxy bulges 
are p robably largely in place at z ~ 2 (jGenzel et al.l 
2008). Accretion onto the black hole at the current rate 
(0.4 M.Edd ~5x 10~ 2 M Q yr _1 ) will bring these galaxies 
onto the local MBH-Mbuige relation (Mbh ~ 10 s M Q ) in 
~ 2 Gyr. Similarly, galaxies in the intermediate mass 
bin will acquire SMBHs with final masses ~ 2 x 10 7 M Q 
in < 1 Gyr. The evolutionary link between these final 
SMBHs at z = (M B h > 2 x 10 7 M ) and their puta- 
tive z ~ 2 host galaxies (M* > 10 10 M ) is supported 
by the similar space densities of these two populations 
(~3x 10 ~ 3 h% n Mpc" 3 and 2 x 1CT 3 MpcT 3 , re- 
specti vely: [McLure fc Dunlod 12004 iDaddi et all l2004bl . 
I2005t iReddv et alj 120051 : iGrazian et all I2007D . More- 
over, the final black holes (Mbh ^ 2x 10 7 M Q ) are 
found at z = i n elliptical and bulge-dominated spi- 
ral g alaxies (e.g. iTremaine et al.1 120021 : iMarconi et al.l 
2004), consistent with the probable descendants of the 
z ~ 2 star-forming galaxy population (|Genel et al.ll2008l : 
IConrov et al.l l2008h . This evidence thus confirm the 
plausibility of the rapid bulge formati on seen at z ~ 2 in 
the SINS galaxies (jGenzel et al.ll2008| ) being followed by 
a few Gyr of rapid SMBH assembly, ultimately resulting 
in spheroids and bulges that obey the local MsH-Mbuigc 
r elation. 

lAlexander et all (|2008l ) found similar delayed SMBH 
formation in the sub-mm galaxy population at z ~ 2 
(see Figure [5]), suggesting that the time lag between 
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black holes and bulges may be a common phenomenon 
in rapidly forming galaxies at high redshift. However, 
this trend is not universal; quasars and radio galax- 
ies at similar redshifts are suspected to lie above lo- 
cal black hole scaling relations, with black holes that 
are over-massive for their host bulge s by up to and ex- 
ceeding an order of magnitude (e.g. IWalter et al.lj2004| ; 
Shields et al.ll2006t iMcLu re et al.l l2006l; iPeng et aLlbooa 



Maiolino et alj 120071 but see also iShields et afl T2003) . 
These systems populate the highest mass end of the black 
hole mass function, with black hole masses of > 10 s — 10 9 
Mq already in place at z ~ 2. It may therefore be that 
black holes in these different mass/activity regimes grow 
in very different circumstances and consequently relate 
to their bulges very differently. If this is the case, the 
challenge is then to locate the mechanism(s) that bring 
these varied high-redshift formation processes together 
into the black hole scaling relations observed at z = 0. 

5. CONCLUSIONS 

In stacked, average spectra of SINS z ~ 2 star-forming 
galaxies, we have detected broad emission underneath 
the much brighter narrow Ha and [Nil] emission lines. 
This broad emission accounts for ~ 30% of the total Ha 
luminosity of these galaxies and can be parameterized 
equally well with a single broad Ha line ( "broad line" of 
FWHM ~ 1500 km s -1 ) and with a two Gaussian fit to 
both the permitted and forbidden lines ("broad wings" 
of FWHM - 550 km s" 1 ). The luminosity and FWHM 
of the broad component increases with increasing galaxy 
mass and therefore with SFR. This broad component is 
found both in known AGN and in stacked spectra of sys- 
tems that have not been previously identified as AGN. 
There is some evidence that the broad emission is more 
luminous in galaxy centers, as opposed to in the outer 
regions, but the significance of these detections are low. 

We cannot empirically determine whether this broad 
emission is due to high-velocity galactic winds and the 
associated shocks or to the BLR emission of AGN. In the 
former case, we find that simple scaling arguments show 
that the luminosity and FWHM of the broad emission 
can plausibly be accounted for via shocking of the am- 
bient interstellar media from supernovae-driven galactic 



winds. These winds would then be ejecting matter from 
the host galaxy at rates slightly exceeding the star forma- 
tion rate, in keeping with expectations from the metal- 
licity evolution of these galaxies and with ultraviolet in- 
terstellar absorption- line studies at similar redshifts. 

On the other hand, the broad emission may be gener- 
ated in a BLR; in this case, we can estimate the black hole 
masses and luminosities necessary to fuel the observed 
emission for each of three galaxy mass bins. We find 
that the measured SMBH masses correlate with the host 
galaxy masses, as expected from local scaling relations, 
but that the SMBHs are significantly under-massive for 
their bulges when compared with local relations. While 
this result has large uncertainties, it is consistent with the 
emerging picture of galaxy assembly at z ~ 2, in which a 
gas- rich disk fragments into large (> 1 kpc) super-star- 
forming clumps that then migrate into the galaxy center 
on Gyr-timescales to form a nascent bulge. The bulge 
would then form first through this process and only later 
completely assemble its black hole. 

The obvious direction for future research is to deter- 
mine the source of the broad Ha emission in high-redshift 
star-forming galaxies. This will most likely require de- 
tailed examination of individual galaxies. Among the 
diagnostics that will be useful for this task are compar- 
isons with X-ray data and deep integrations in rest- frame 
UV/optical wavebands to spatially resolve e.g. UV inter- 
stellar absorption lines and broad Balmer emission. 
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